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ABSTRACT 


WV  ro:i^i(ier  tlie  pruhleni  of  diri'ction  tiiufiiis  (DF)  of  inuhiplc.  cofrequcncx' 
narrowhaiul  sisnals  with  a  phased  aiiteiiiia  array  when  direcrloii-de])endcnt  arra\- 
calihiratioii  errors  as  well  as  thermal  noise  are  present.  Lower  hotintls  on  tlie  varianee 
at  mihiased  estimators  for  DF  are  derived  under  two  difff  rtmt  t\  pes  of  si”;nal  model'; 
completelx'  unknown  (i.e..  gtmeric)  signals  and  unknown  constant-envcloiie  sijtnais. 
In  both  models,  the  romple.x  amplitudes  of  the  sigtials  a'e  modeled  as  unknown  pa¬ 
rameters.  W'e  deri\  e  and  evaluate  lower  bounds  oii  DF  of  these  signals  when  romjtle.x 
Ciatissian  arrtiv  errors  and  thermttl  noise  are  present.  In  our  nutnerical  e.xamples. 
the  Itound  for  generic  signals  temled  to  decrease  with  increasing  interference  power 
aiifi,  for  closely  spaced  signals,  became  more  optitr.istic  when  suiall  airay  errors  were 
iidded.  \\'ith  a  suflicicntly  large  signa!-of-interest  (SOI)  array  signal-to-noise  jtowi 
ratio  (.ASNR)  and  tmmber  of  looks,  the  bound  numerically  approached  the  bouini 
on  FtF  of  multiple  generic  signals  with  array  errors  but  no  thermal  noise  (i.e..  the 
multiple-generic-signal.  arra\-errors-onl>  bound).  Tlu'  hitter  bound  depended  (ui 
the  signal  separation  and  tlie  power  of  tlie  array  errors.  It  wa.s  indepci.ucnt  of  die 
signal  wai'eforuis  and  powers  ai.d  was  largest  at  moderate  signal  sejiarations  t'.g.. 
0.:i  lieamwidths  .  The  latter  bound  approached  the  single-signal  array-errors-only 
bound  as  the  signal  separation  approached  zero. 

The-  results  for  consttint-envelope  signals  showed  tlnit  the  bound  with  arrav 
errors  and  thermal  tioise  was  more  optimistic  than  the  analogous  genedc-signal 
bound  anrl  had  little  dependence  on  signal  separation  and  interference  power.  Thi' 
w;is  tilso  the  case  in  the  absence  of  direction-dependent  errors.  The  imprcnemetits  in 
tile  bound  were  more  significant  at  small  arra\  errors,  low  .ASNR.  small  numbers  of 
looks,  and  moderate  signal  seitaral ions.  .-Xt  all  signal  separations,  with  a  sufficientli 
liiali  SOI  .4S.\R  and  number  of  looks,  the  bound  niimericai’ 
signal,  otie-look.  arr;iy-errors-oul\  bound  (the  latter  bouiu 
and  constant -emelope  signalsi. 
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PREF/  CE 


lo  ix'rform  direction  finding  (DF)  of  imiltijdo,  cofretjuency  narrowltand  signals 
with  a  phased  anttnina  array,  it  is  necessarv  to  know  the  arra\ -response  pattern  as 
•I  function  ot  tfie  signal  direction.  In  practice,  tftis  function  is  known  only  approx- 
imatel\-.  1'  is  tlierefore  appropriate  to  model  the  array-response  paitern  vector  as 
a  perturbed  versiftn  of  a  deterministic  function  of  tiie  signal  direction  (i.e..  the  as¬ 
sumed  array-response  \-ector).  It  is  n.seful  to  decompose  th('  perturltations  into  two 
components,  a  direction-independent  component  and  a  direction-dependent  compo¬ 
nent.  The  direction-independent  component  is  commonly-  modeled  as  multiplicative 
errors  that  are  the  same  for  all  of  the  signals  in  each  recci\er  channel  le.g..  in¬ 
ternal  receiver-calibration  errors).  The  direction-dependent  component  is  due  to 
errors  that  are  different  for  each  of  the  signals  in  each  receiver  channel  le.g..  ei- 
rors  due  to  array-respon.se  pattern  calibration  residuals  and  near-field  mnhipati) 
The  direction-dependent  component  is  modeled  as  adding  errors  to  the  assumed 
array-response  vectors  of  the  individual  signals.  In  this  report,  we  shall  limit  uir 
consideration  to  the  directiori-dependetit  componer,t  of  the  errors. 

Direction-dependent  array  errors  can  be  an  important  limiting  factor  in  ilie 
performance  of  DF  algorithms.  To  help  determine  whether  current  algorithms  are 
performing  near  the  intrinsic  limitations  of  the  problem,  it  is  useftil  lo  compute 
lower  bounds  on  the  variance  of  unbiased  DF  estimates  in  tlie  presence  of  direction 
dependent  e  rors. 

in  Section  1  of  thio  report,  we  derive  lower  bounds  on  the  variance  of  unbiased 
estimators  for  DF  of  unknown  generic  and  constant-envelope  signals  when  inde¬ 
pendent.  complex  Gaussian,  direction-dependent  array-calibration  errors  as  well  as 
complex  Gaussian  thermal  noise  are  present.  (Note  that,  in  a  real  er.vironment.  the 
arrav  errors  might  be  correlated  ant:  nna-element-to-antenna-element  and/or  signal- 
direction-tosignal-direction.  which  could  impact  the  performance  achievable.)  Nu¬ 
merical  results  for  some  cases  of  interest  are  presented  in  Section  2.  Conclusions 
and  directions  for  ♦'uture  work  are  presented  in  Section  3  and  Section  4.  respectively. 


1.  SIGNAL  MODELS  AND  LOWER  BOUND  DERIVATIONS 


1.1  NOTATION  AND  SOME  MATHEMATICAL  NOTIONS 

I^eal  and  complex  scalars  are  denoted  by  lowercase  Roman  or  Greek  letters.  \'ectors  are 
(ienoted  h\  lowercase  Roman  or  Greek  letters  with  an  arrow  on  top  (e.g..  /).  Matrices  are  denoted 
li>  ui)percas('  Roman  or  Greek  letters.  Exceptions  are  the  symbols  used  to  denote  the  dimensions  of 
vectors  and  matrices  as  well  as  the  symbols  used  to  denote  the  number  signals,  antenna  elements. 
oh.ser\at ion.s  (i.(\.  looks),  and  unknown  parameters,  .lamely  S.  M.  .V.  and  P.  resirectivcly.  The 
letters  .s.  in.  n.  and  p  are  used  to  index  over  tlu'  signals,  antenna  elements,  looks,  and  unknown 
p;ir;nnt'ters.  respectively.  The  dh  column  of  the  matrix  A'  will  be  denoted  by  A'.,,  the  dh  row  of 
thi“  matrix  .\  (  written  as  ti  column  vector)  will  be  denoted  In'  A,.,  and  the  ijth  component  of  the 
matrix  A  will  be  denoted  by  iA'),j.  The  dh  component  of  the  vector  T  will  be  denoterl  by  (Ji,. 
l  ilt'  re;il  and  imaginarv  pjirts  of  a  quantity  A  will  be  denoted  by  'R(A)  and  T(A').  respectively. 
I  he  trtuispose,  complex  conjugate,  and  Hermitian  (i.e..  complex  conjtigate  transpos(')  of  a  matrix 
.\  are  denoted  b\'  A’ ^ ,  A',  and  A^^  respectively. 

We  define  f,  as  ttie  vector  whose  ith  component  equals  1  and  whose  remaining  components 
ecpial  ti,  1  V  as  the  .V-vector  whose  components  all  equal  1,  /y  as  the  .V  A  identity  matrix.  v 
as  the  M  ■  .\  imitrix  whose  components  all  equal  0.  and  E,j  as  the  matrix  whose  tjth  component 
ecjuals  1  <ind  whose  remaining  components  equal  0.  We  shall  need  the  identities 

I.T  i  ABi  -  tr  { BA}  f  ] .] ) 

2.  r.r|,  1  .  111  .  and 


tr  ( .4.,ej  ) 

(12.  Eq.  1.28] ) 

r.r(cjA,) 

(Eq.  1.1) 

tr  i(.4)j,; 

(1.2 

1.1.1  The  Kronecker  Product 

I  he  Kronecker  product  of  the  I  '  ,J  matrix  A  and  the  K  x  L  matrix  B  is  defined  as  the  IK  >  ,J L 
matrix  2.  Sec.  2.2 

(A)uB  {A)i2B  .  •  {A)uB 
iA)'2lB  {A)22B  ■■■  {A)2jB 

{A)iiB  {A)i2B  (A)fjB 


It  has  the  properties 


A  (aB)  =  (a.4)  B  =  q(.4  B) 


(1.4) 


(,‘i.  Eq.  2..')  and  a  trivial  extension), 

(.4  Bf  =  /?" 

(.1,  p.  B]. 

(.4  ■;  Br^  =  ,4-‘  B-‘ 

(  2.  Eq.  2.12  ).  if  the  invenses  exist,  and 

(  .4  B)(C  D)  =  .dC  BD 

(  2,  Eq.  2.11  ).  where  C  is  a  J  x  M  matrix  and  D  is  an  L  x  A'  matrix. 


(1.5) 


AAV) 


1.1.2  The  Hadamard  Product 

The  Hadamard  product  of  the  I  y  J  matrices  A  and  B  i.s  defined  a,s  the  /  x  J  matrix 


AAB  1^' 


(.4)ii(5)ii  [A)i2(B)i2 

{A)2\(B)2\  (.4)22(5)22 


{AhAB)u 

(A)2j[B)2J 


{A)i\{B)i\  (A)i2(B)i2  {A)ij{B)ij 
The  Hadamard  product  is  obviously  commutative  and  associative,  so  that 


and 


AOB = BOA 


{A0B]0('  ==  AOIBOC). 


W'v  shah  fteed  the  identities 


f.4i;/l)'  =  a’ OB’ 


{ah’  ).-.{c(P  )  =  {dOr)U>^d)‘ 


(1.8) 


(1.9) 

(1.10) 

(1.11) 

(1.12) 
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and 


(E^rA):n(BE,t)  =  (12.  Eqs.  1.27  and  1.28:) 

=  (Eq.  1.12) 

=  {A)rt(Bpi-pJ 

=  (-4)r((/^)g.,Eg(  ([2.  Eq.  l..a;).  (El.'i) 

1.1.3  The  Khatri-Rao  Product 

The  Kliatri-Ilan  product  of  the  M  x  :V  matrix  .4  and  the  P  y  N  matrix  B  is  defined  by  the 
MB  ■  A'  matrix 

.4  .  Z?  '=  [  .4  1  5,,  .4.0  T  B.o  •  •  •  -4.,v  £  B  ]  .  (1.14) 

\\V  shall  need  the  identities 

(.4  :■  B)(C  I  D)  -  {AC  7  BD)  (jl,  Eq.  1.5.9!)  (11^) 

and 

(.4  '•  Bf{C  I  D)  =  (A^ C)n(B" D).  (1.16) 

The  latter  is  proved  by  noting  that  the  ijth  element  of  the  left-hand  side  is  equal  to 

{A.,P  B.,f(C.p-:  D.j)  =  (A^  Bp(C,  &  Dp  (Eq.  1.5) 

=  iA^,CPrZ{B^Dp  (Eq.  1.7) 

=  {A^ChjCPB^np 
=  (.4"C)q(B''D),,. 

We  shall  also  need  the  identity 

(.4  •  BpiC  D)[E  ■  F)  =  {A  -  BfiCE  •  DF)  (Eq.  1.15) 

=  A^CEOB^^DF  (Eq.  1.16). 


(1.17) 


(1.18) 
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1.1.4  Matrix-Inversion  Identities 


W'e  shall  need  the  matrix  identity 

:.-l  fiC  'Zlj  ‘  =  .d'*  -  DA'^B'r^DA-^  (1.19) 

(assuming  the  necessary  inverses  exist)  a.  Eq.  AM],  and  the  following  result  on  the  inverse  of  a 
partitioned  matrix  a,  E(:.  .‘V.'i,.  Let  .4  be  an  invertible  square  matrix  partitioned  as 

t'l- 

*^21  -H22 

with  .-l.)_;  invertible.  Then  the  upper-left  block  of  .1  .-1'',,.  is  givtm  by 

•■5  'll  (*^1!  -  -dlJ.doo  (1.21) 

(the  inverse  indicated  on  the  right  side  of  the  equation  can  be  shown  to  exist  under  tlie  given 
conditions),  and 


.-r' 


A- 


11 


— -doo  Any  A 


-) 


11 


-t-  .doo  '.42l.d'^  J,.di2.d2.r' 


(1.22) 


Note  the  nolational  distinction  between  the  block  of  .d  ^  .d  and  the  inverse  of  the 
block  of  .d. 


1.2  DATA  MODEL 

1.2.1  Array-Response  Model 

The  signal  environment  will  be  modeled  as  S  narrowband  cofreqnency  signals.  The  complex 
amplitude  of  the  .sth  signal  during  the  nth  observation  (i.e.,  look)  will  be  denoted  by  c,,.,.  The  c„,, 
will  represent  an  unknown  generic  or  constant-envelope  signal.  The  details  of  these  signal  models 
will  be  described  in  Stibsections  1.2.2  and  1.2.3.  We  define  the  5-vector  of  complex  amplitudes 
durin<£  the  nth  look.  c„,  by 


<''n  ='  [  fV,i  Cn2  ■■■  CnS  J  .  (1.23) 

the  concatenation  of  these  vectors,  c.  by 

d'*-  ^(c„  Xc;,),  (1.2d) 

n=l 


4 


and  the  A  -  .S'  niatri.x  of  ruinplex  amplitudes,  (’.  liy 


Cll 

(‘l-j  - 

C 1  s 

Cil 

<''2S 

c.vi 

C\  2  ■ 

(‘VS 

We  will  assume  that  .S'  and  that  (’  is  of  full  rank.  .S'. 

The  signals  are  ubservetl  using  an  arra\‘  of  M  antenna  elements.  'I  ht'  ideal  (i.e..  errur-lree; 
arra\-response  \'ector  uf  tlu'  .sth  signal  will  he  dmioieil  by  m.  and  the  errnr  in  this  \ect(.)r  will  lx* 
denoted  by  Define  the  romjjlex  M  •  .S  matrices  ol  ideal  iirray-rt'sjionse  vectors  and  their  (.‘rrois. 
tciid  E.  I'.y 


-  def  r 
-  :  Cl 


and 


■  c;  E,  ...  E. 

—  1  si  si  S,'- 


The  resjronst'  of  the  arrtiy  during  the  /ith  look'.  E„.  will  be  modeled  as 

-  liyf  ,■  I  •  3  1  r  - 

-71  —  (1  —  IC;)  S  . 


( i.2(;i 


(1.27 


(1.2  m 


where  k„  denotes  the  thermal-noise  vector  corrupting  the  observation.  The  vectors  {Tv}  ttnd  {E,, ’ 
will  be  modeled  as  seciuences  of  independent,  zero-mean  complex  (circular)  Gaussian"  random 
vectors  with  covariance  matrices  f’/y;  tind  (T"/y/.  respectively.  The  {Tv}  <md  {E,,}  are  assumed  to 
be  statisticallv  independent  of  each  other. 

Define  the  concatenated  thermal-noise  vector.  E.  by 

.\ 

E  ’=  ^(E„  ■  H„)  (1.2'n 

n-  I 

and  the  concatenated  obser\'ation  vector.  5.  by 

n  -  1 

'  .A  look  (or  observation)  is  defined  as  a  sinndtaneous  vector  snapshot  of  the  baseband  complex 
am[)litudes  at  the  outputs  of  the  M  antenna  elements  composing  the  array. 

~  In  the  rejtort.  we  shall  often  use  the  terminology  'complex  Gaussian'  in  phice  of  ’complex 
circular  Gaussian'  for  sini()licity. 

7) 


.\ 


n  1 

A'  >'  ,\ 


I 


n  \ 


S  .V 


5Z5Z  ^  ^  ->v  *  ^  (^‘^1- 


i::r..:,(r*E)..: 


,-■ '  1 

C  •  (\' -  H)  U  -  K  (Kq.  1.14). 

^hu^  the  veetd!'  5  is  a  ooiuplex  Gaussian  rancioni  .\/A-vector  with  mean 
/7.=  (C  ■  l')l's 

and  covariance 

.\  E{'{c  ■  E)u-  -  K  ;ri'(r  •  E]“  ^  h'^]} 


E{Y,iC.,  :•:  ^J}^E{Rh>^) 

•s  -  1  > '  i 

s'  v' 

'•  ^  E{RR^')  (Eq,  1,- 


Is'- 1 

S'  s' 


s=  Is'-  1 
5 

=  ;•- f-/A,i-tT‘(/v  eIm) 


,<=1 


=  {a‘E\  -rCC^<]  y.  Em  (  2.  p.  16,). 


Define  .\  \'  by 


then 


.\.v  =  rr'^/.v  -  c'^GC": 


(l.:5(b 


(i.;5i 


(1.32) 


(1.33) 


(1.34) 


A  -  A.v:-. /.u. 

W'e  will  assume  that  either  rr-  >  0  or  that  both  :V  =  5  and  e‘  >  0  so  that  Ay  is  of  full  rank.  .V. 
Then,  using  Eciuation  (1.6). 


.\  '  =  Im. 


(1.3.3) 
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(fT'/v  *  f-CC’’)-' 


with  Ay'  heing  givyn  b>’ 


'^{/v  ^  — rr")  - ' 

rT“ 

Ai/v-  -.cc"y' 

O'-  Cf'- 

1  f~  “ 

—  7v  -  (Eq,  1.19) 

(7-  (T-  rr- 


(l.iUI) 


whon  (T“  ■  0,  a;  d  h\' 


Ay'  -  il.;57) 

f  - 

when  rr“  =  0,  A’  =  5.  and  f“  ,>  0. 

The  id»'al  array-resfjonse  vector  of  the  sth  signal,  t^.  i'^  assumed  to  be  a  differentiable  function 
of  the  one  dimensional  direction  of  the  sth  signal,  which  we  denote  by  a.,.  (The  a.ssumption  of 
a  one-dimensional  space  of  signal  directions  is  made  only  for  simplicity;  the  method  of  analy.sis 
could  be  extended  to  a  multidimensional  space  of  signal  directions.)  The  real  5-vector  of  signal 
directions,  u.  is  defined  bv 


-  def  r  ^ 

a  -  ,  u-:  ■  •  •  a-5 


(1.3») 


We  denote  the  derivative  of  c,  with  respect  to  a*  (evaluated  at  the  true  value  of  u.,)  by  r.,  and 
define  the  matrix  V'  bv 


def  r  u  U  U 

\  =  r'l  V2  •••  I’.s 


(1.39) 


1.2.2  Generic-Signal  Model 

In  the  generic-signal  model,  the  complex  amplitudes  of  the  signals  are  modeled  as  being  un- 
knowii.  More  specifically,  c„s  will  be  modeled  as 

Cn.s  =  (1.40) 

where  and  Ons  ate  the  unknown  magnitude  and  phase,  respectively,  of  the  ,sth  signal  during 
the  nth  look. 

We  define  the  real  5'-vector  of  signal  magnitudes  during  the  nth  look,  d^.  by 


—  def  ] 

(iji  —  (lyjl  ‘  ’  ^nS 


the  roncalcnation  of  these  vectors,  a.  by 

(1.42) 

ii  =  i 


the  real  vector  of  pha.ses  during  the  nt)i  look.  by 


r  dcf  ' 

Oji  Ojil  Oji2  '  '  ■  OfiS 


(1.43) 


and  the  coucateuation  of  the.se  vectors,  o.  by 

(1.44) 

We  also  define  the  complex  5-vector,  u'n-  of  phases  (in  exponential  notation)  during  the  tith  look 
V)\- 


(i-'rj  — 


(1.45) 


and  the  concatenation  of  these  vectors,  y,  by 


.V 


n=  I 


The  vector  of  unknowns,  f.  is  defined  by 


_  def 
r  = 


-T 


(1.46) 


(1.47) 


1.2.3  Unknown  Constant-Envelope  Signal  Model 

In  the  unknown  constant-envelope  signal  model,  the  complex  amplitudes  of  the  signals  are 
modeled  as  being  of  unknown,  but  constant,  magnitude  and  of  unknown  phase.  More  specifically. 
c„,,  will  be  modeled  as 

c.„.;  =  (1.48) 

where  a,,  is  the  unknown  magnitude  of  the  sth  signal  and  (/)„.,  is  the  unknown!  phase  of  the  sth 
signal  during  the  nth  look. 

We  define  the  real  5- vector  of  signal  magnitudes,  q,  by 


8 


_  lief  " 

( 1  -  I  a  1  ( i  j 


ami  0,1  aiul  a>  in  Sectiun  1.2.2.  1  Im  von  or  of  imkiKJwiis.  r.  is  dofineci  !)>■ 


...Icf 

r 


,7'/ 


’  7 


(i-ani 


1.3  The  Cramer-Rao  Bound  on  DF 
1.3.1  Tlie  Cramer-Rao  Bound 

Define  as  tlu'  row  \eetor  whose  /th  eoniponem  is  equal  to  ■  Let  5  he  a  ramloin  \'ector 
wliose  [jo.ihahilit y  den.sitx'  funrtioti.  />{::  r).  (lepends  upon  an  unknown  Z.-diinensional  real  jraraiuelei 
\eetur  r  with  true  vtdue  />.  LtU  r  dt'uote  ;ui  estimator  of  r  (F  is  ti  function  of  5)  and  define 


-  ,  -  ,  dl  l'  ... 

nitr/]  -  Lr.  (/•). 


wliere  £V.  denotes  expecttttion  witli  resitect  to  the  proltahility  densitv  function  /t{5;  r,  )■  The  covari- 
;uic('  of  the  estinuttor  r  when  r  =  F/  is  defined  hy 

I:):.  {  F-  m ( F, )'  F  -  iTi{  Ft  f ( 1 , .') 2  i 

We  sa\’  that  r  is  different ialh-  unbiased  for  (F)/  at  D  if 


(  Urn  I  F}  \ 


(LTD 


and  tfi;tt  F  is  differentially  unbiased  ttt  F,  if 


dw{  F) 
_  ----- 


-  //,. 


(l.o-D 


Ilefiiie  tin'  Fisher  information  matri.x  by 


F  E 


dlnp(.;';F) 

()F 


din p{F:  F] 
dF 


The  Crtimer-Hao  bound  states  th;tt 


H 


(l.ah) 


(l.'xi) 


!> 


where-  >  nieans  that  the  ciift’ercncc  between  the  matrices  is  positive  semidefinite  !3.  Th.  2.7.3’. 
If  F  is  differentially  unbiased  for  (r)/  at  ry.  then 


var  {r]i[ 


iR)u 

ej 

Rfi 

( dll 

\( 

dll) 

[1 

df 

rf 

F-\- 

[F 

,  ,  dm 

^  '  Hr 


e/ 


'■'/  i  J 


(i2.  Eqs.  1.16  and  1.17  ) 


l..')7: 


In  this  case,  the  bound  on  the  \ariance  of  the  /th  unknown  parameter  is  gi\en  simply  In  the  /th 
diagonal  element  of  f  ”*.  In  the  sequel,  for  any  parameter  of  interest,  (r)/.  we  shall  only  consider 
estimators  which  are  differentially  imlhased  for  {i^i  at  E.  For  an  estimator  which  is  differentiall\ 
tinbiased  at  o.  Equation  (1.56)  reduces  to 


R>  f-'. 


(1.5M 


Lastly,  we  note  that  locally  unbiased  e.stimators  of  (f)/  at  fi  (i.e..  estimators  for  which  'm(r)  / 
=  (F)/  for  all  Fin  some  neighborhood  of  F, )  are  always  also  differentially  unbiased  estimators  of  (F); 
at  F.  This  is  shown  by  noting  that 


/  Omif) 

y  di  f-f/ 


h 


Of 


I 


ih),. 


</• 


(1.5!)! 


1.3.2  Generic-Signal  Case 

I’sing  Erjuation  (1.47).  the  symmetric  matrix  F  may  he  partitioned  as 


*-  uu 

F- 

^  uo 

I/O 

/•' 

— 

F- 

^  a  u 

F-- 

a<i 

r 

ao 

.  ^0,7 

F^. 

Oo 

F- 

oo 

10 


^  U/i 

^  ud\ 

‘  lias 

F^- 

iiOi 

-■  F-- 

U0\ 

^  li  \  ii 

^  <1  j  1 

F-  -  • 

a,o, 

•  •  / 
aj  O.v 

F-  - 

.  .  ^ 

■*  o  .va  v 

F.  -  ■ 

a  vOi 

^Q.vO.v 

F-  -- 

Oi  O 

^  • 
0 1 « 1 

-  •  F-  - 

OiO.V 

OlOl 

. .  f.. 

OlO.V 

F-  .  F-  -  ■  F-  .  F-  - 

L.‘_v  ti  C'.v«i  o.va.v  c  .vOi 


where  F is  detined  b>- 


F  -  r 

O  .V  .V 


(].()()) 


f  h:  =' 


d]np{::r)\  dlnp{z:/^ 


c)3 


03 


1.61) 


■Since  F  is  syninietric.  we  need  only  explicitly  compute  the  diagonal  and  lower-left  subljlocks  of  F. 

The  ol)scrvation  is  a  complex  Gaussian  random  vector  with  mean  /7  and  covariance  .\.  For  this 
case,  it  follows  from  1,  Eqs  2.6  and  2.23^  that  the  /jth  element  of  F  is  given  by 


F,j  =21? 


0^ 

Or, 


Oi-J 


-  tr 


Or, '  Oi-j  ) 


(1.62) 


where  .\  and  the  derivatives  are  evaluated  at  the  true  value  of  the  unknown  parameter  vector. 
Explicit  expressions  for  the  submatrices  of  F  are  given  in  Appendix  A. 


1.3.3  Constant-Envelope  Signal  Case 


Using  Ffquation  (l.-aO).  the  symmetric  Eisher  information  matrix  F  may  be  partitioned  as 


^  u  (7 

^  liO 

^  Oli 

^  ad 

^0-0- 

F-. 

F-, 

F-r 

Ou 

Oa 

oo 

^  no 

^  t7d 

/■.- 

UO] 

•  ■ 

uO\ 

^  d  Ii 

F-  - 

^  (»f» 

f\^  ■ 

■  ■  'f'-g 

aO.v 

F-  - 

0\  II 

F  -  - 
0 1  <> 

F-  -  ■ 

OiC'i 

p 

OiO.v 

(1.63i 


r  -  -  r  -  .  r  ■  ■  •  t  ~  ~ 

L  0\u  Cfv'i  O.vCm  O.vCfv  J 

.As  in  the  generic-signal  case,  tlie  observation  is  a  complex  circular  Gaussian  random  vector  with 
mean  p  and  covariance  .\.  so  that  the  ijth  element  of  F  is  again  given  by  Equation  (1.62).  It  is 
shown  in  .Appendix  B  that  the  information  matrix  for  the  unknown  constant-envelope  signal  case 
is  easily  constructed  from  the  information  matrix  for  the  generic-signal  case. 
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1.4  A  BOUND  ON  DF  OF  STOCHASTIC  SIGNALS 


The  bounds  on  estimation  of  u\  derived  in  Section  I.. 1.2  and  Section  l.l.l  depend,  in  general, 
on  the  specihe  complex  amplittides  of  the  signals  (as  well  as  their  directions)  and  apply  to  estimators 
which  aie  differentially  unbia.sed  for  U]  (see  Section  1.1.1).  In  this  section,  we  will  use  this  bound  to 
obtain  a  lower  bottnd  on  tlie  \'ariance  of  a  more  restricted  class  of  estimators  for  the  case  where  the 
complex  amplittides  of  the  signals  have  a  probabilistic  distribution.  This  bound  has  the  desirable 
property  of  not  depending  on  the  specific  complex  amplitudes  of  the  signals. 

Partition  the  parameter  \ector  r  into 


W'e  say  tluit  d is  gloltally  different ially  unbias(>d  for  with  respect  to  r'  and  different iall\'  unbiased 
with  respect  to  P  at  if 

f  (bn  I  Pi 


holds  for  ;dl  possifih-  ir').  Note  tluit  this  im[)iies  that  ( tn( )  i  is  constant  with  n'spect  to  r". 
For  th('  remainder  of  this  section,  we  will  refer  to  this  (juantity  as  mi(f^). 

In  this  section,  we  consider  estiniiitors  of  i/i  which  are  globtilly  different iall,\  unbiased  with 
respect  to  the  signal  complex  amplitudes  as  well  as  differentially  unbiased  with  respect  to  ilie 
signal  directions  at  .Note  that  the  reciuirement  of  global  differential  unbiasedness  with  respect 
to  t'le  complex  amplitudes  essentially  metins  tluit  the  estimator  still  treiits  the  complex  annilitude' 
as  being  unknown  parameters  rather  than  ;is  having  some  known  jirobabilily  distribution.^ 

I.et  P'  denote  the  parameters  associated  with  the  com]rlex  amplitudes.  These'  parameters  are 
noM-  assumed  to  ha\e  some  fixed  probaliility  distribution.  Let  t,  denote  the  remainder  of  the  random 
\ariaf)les  in  the  problem.  Let  /Pc’.P')  denote  the  joint  jirobtibility  density  fitnction  of  C  and  r"'. 

'  .A.!!  estimator  which  'knows'  the  probability  (listribution  of  the  complex  amplitudes  is  said  to 
be  differetitially  unbiased  for  t/i  at  r,  if  it  satisfies  Equation  (l.a.'I)  with  /  =  1  and  p  =  /?/  (i.e..  the 
unknown  jiarameter  vi'ctor  does  not  include  the  complex  amplitudes),  where  the  expectation  in  the 
dehnition  of  w  Ecpiation  ( 1 .51 )  is  taken  with  respect  to  the  distribution  of  the  complex  amplitude's 
as  well  as  the  other  random  f|uantities  in  the  problem.  It  can  be  shown  that  any  estimator  of  ui 
which  is  globallv  different iall>-  unbiased  with  respect  to  the  signal  complex  amplitudes  as  well  as 
differentially  unbia.sed  with  respect  to  the  signal  directions  also  satisfies  this  condition.  Howevi'r. 
the  cotivers('  statement  is  not  true. 


I.ti5,: 
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Let  y  he  an  unlhased  (in  the  st-nse  given  above)  estimator  of  the  direction  u\.  Tlie  mean  of  y  vvitli 
resiteot  to  tlie  joint  distril)nti(jn  of  C  mid  is  given  by 


ypiC-  )  dC 


■IQ  l.iy 


yp[C,  f")  dC 


»i\{r,)p(r")  dr' 


(1.6G) 


wlu-re  Z  denotes  the  domain  of  d-  denotes  ihe  domain  of  r".  and  the  second  equalit\'  follows  by 
the  comment  after  L'qtir'ttion  (  I.rioj.  Let  //'  denote  the  variance  of  y  witli  respect  to  p((,.r').  Let 
ss  denote  tin'  ( 'ranier-R;\o  liound  (obtained  in  Section  1.15.2  or  Section  1.15.15)  on  the  estimation  of 
u\  tor  tile  parameter  vector  r".  W'e  have 


I  11  >11 1 1 1^1  '  "/'( s .  f*'  i  dd  dZ' 

j  I  y  HI  1 1  /■' )  ~pi  s'  r"  i  ds'  •  />(  r' )  dr 

V/  .  j  41 


p>  r  I  dr 


(Ltd 


1  hi'  i'  the  de'ired  lnwer  bniiiid  on  the  variance  ot  (‘st imtUttrs  (.li  (/[  which  are  globtilly  diH'erentiallv 
imfiiased  with  rt-sfierf  Ui  the  >iyn<'il  roir.plt's  at/iplil tides  ant]  flitferf'itliallN  nttbiased  witli  respect  to 
the  signal  directions.  In  the  se(|nel.  we  will  evalintte  this  itonnd.  which  we  term  tlu'  'integrtited' 
bniind.  We  wi!i  make  a  few  e.xcepfintis  to  this.  Itectnise.  in  st.ime  cases,  the  Cramer-Rao  bounds 
derixed  in  St'ciioii  1.15.2  and  Section  1.15.15  turn  out  to  be  numerically  itidepetuient  of  the  specific 
complex  amplitude.-  of  the  signal'  Iti  such  cases,  tlie  integrated  bound  in  Kquatioti  (l.h7i  is 
nnmericalls  identical  to  the  (  ramer-Rao  bounil  which  is  obtaiiu'd  for  an\'  itllowable  set  of  comple.x 
aiiipiit  ndc'.  In  'iich  case',  we  find  it  useful  to  contiinn'  to  reh'r  to  the  rele\’anl  Cranier-Hao  bound 
rather  than  the  integrtited  bound.  In  wh.it  lollows.  the  term  Crauier-Rao  bound  will  be  used  wlu'U 
and  onl\  wlien  i  he  integrtited  boumi  gi\en  by  lapiatioii  (l  .tiil  is  not  intetided. 


2.  NUMERICAL  EVALUATION  OF  THE  BOUNDS 


2.1  METHODS 


2.1.1  Array  Signal-toNoise  Power  Ratio  (ASNR) 

111  what  follows,  it  is-  useful  to  define  the  .ASNR  of  a  given  sourcf'.  In  our  hound  evaluations, 
we  will  con.dder  onl\  stationary  stochastic  si^iial  sources.  We  define 

l\  L'(  -‘i;  yn  -  1 . ==  1 . .S’).  ('-■!  ! 


v.here  I.  (ieisites  unit lu'nnit iral  expect ;it ion. 

1  .''ing  Lfjuation  (  1.2^)  in  tiie  ahseiu'e  of  ;irra\'  errors  and  interferers.  if  one  .vere  to  a)ipl>'  a 
tdinph’x  .'\/-\'e('to'-  of  weiglits.  i.\  to  the  tuitenna  element  ontituts  at  the  tuh  look,  one  would  olitain 
an  arra;.  output  of 


The  first  term  in  L(|uation  (2.2)  can  be  identified  a.s  the  signal  component  of  the  array  output  aiid 
the  second  as  the  noise  component.  The  signahto-noise  power  ratio  (SNR)  at  the  array  output  can 
t  hus  he  defined  a.s 


SNR,>, 


{jcf  ill  tC  t  '  ** ) 

E(or^/K„r-) 

_  :  ir  t  V  ’  'EjlCns'^) 

f)s  itr^'rd' 


fT-  IP/'I 


(2.;5) 


H\  the  Seliwar/  ineriuality,  tlu'  weights  which  maximize  the  output  SNR  are  of  the  form  Jvy.  where 
i  In  an  arl)itr;ir\  nonzero  complex  .scalar.  The  resulting  maximized  output  SNR  is  termed  the 


AS.NR: 


,.\SNR  max  SNR.,,,, 

r’ 


(2.4) 


v.iiere,  fipi  ( I ao'enienee,  we  will  assume  in  the  seqtiel  that  the  ideal  a  ray-response  vectors,  the  {ty}. 
are  oi  unit  norm  l  lius  the  .ASNR  of  ;i  given  signal  is  defined  as  the  maximum  SNR  ohtainahk'  for 
th.it  signal  at  the  output  of  the  arrtiv  in  the  absence  of  array  errors  and  interferers. 


2.1.2  Array-Respoii.se  Vector  Model 

rile  ideal  arra\  rt’sjioiise  \v;us  iiiodeled  as  tliat  of  a  imiforin  linear  array  of  ir.otiojiic 

antenna  el('in('nts.  lor  this  array,  the  inth  coinpon'Mil  of  v.,  is  ecpial  to 

>  A  , . ,2.',) 

\  .U 


(i,  p.  72  uIkt* 


Is  the  phase  center  iif  the  array,  A  is  the  wavelenj^th  of  the  carrier,  tl  is  the  inti'releinent  spiiciny,  ol 
the  arra\  (here  (I.2A ).  and  fK  is  the  aneic  of  arrival  of  t  he  .st  h  signal  relative  to  t  lu'  broadside  ))lane 
deiined  as  the  plane  that  jiasses  throiif;h  the  plia.s<*  <'<>nler  of  the  array  and  is  perpt'ndicnhir  to  the 
arra\  axis  (see  figure  2-1  )  .  fau  h  r'^  is  ol  unit  norm,  as  desired. 


f  'j^'iire  2-f  Ihtifiinii  /inear  array  grouit  .'yv. 

We  will  use  beannvidth  (HW  )  units  for  angle  of  arrival,  where  the  angle  <if  tirrival  (d  the  .sih 
siynal  in  BW  is  deiined  by 


,i,.| 


lb 
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The  use  of  B\\  units  is  luotivated  as  follows.  In  Section  2.1.1.  we  detennined  that,  in  the  absence 
of  array  errors  and  interferers.  the  SNR  at  the  output  of  the  array  could  be  inaxiniized  by  selecting 
a  weighting  vector  for  the  array  element  outputs  which  is  a  scalar  multiple  of  the  array-response 
\ector  of  the  source,  i.e.. 


w  =  .dds- 


(2.M 


for  some  nonzero  cumple.x  scalar  .i.  l.'.se  of  this  array  weighting  vector  is  commonly  termed  conven¬ 
tional  beamforming.  If  these  weights  were  held  constant  while  the  source  were  moved  awa\'  from 
direction  to  direction  u.  then  the  output  SNR  of  the  array  would  change  from  the  ASNR.  p.,  a~. 
to  see  Equation  (2.3) 


SNR, 


v(  a ) 

O 

(T“  .:ir  “ 


a- 


(2.!)) 


W'e  thus  obtain  a  'gain  pattern"  for  the  array,  as  a  fnnction  of  the  source  direction.  //.  For  the  linear 
array  which  we  have  defined,  this  pattern  is  given  by 


SNR„u,(u)  = 


p,.  /  sin'rrlw  -  u^} 

CT-  \  M  sin:7r( ;/  -  n.)  M 


(2.1(11 


f).  p.  .  which,  as  expected,  has  its  peak  value  of  p., /rr*  at  u  =  (p.,..  and  also  has  its  first  nulls  (i.e.. 
directions  of  zero  gain)  for  signals  wliicii  are  one  BW  away  from  (i.e..  for  which  u—  ?.■,«  =  2 ).  I'se 
of  B\\'  units  facilitates  performance  comparisons  for  arrays  having  different  interelement  spacings 
e.g..  note  that  Efiuation  (2.10)  is  a  function  of  (/.  but  not  of  (ci  A)  . 

-Array-error  powers  are  given  in  dB  relative  to  the  elements  of  the  ideal  array-response  vectors. 
Since  the  magnitude  of  each  clement  of  the  r.,  is  1;  v' A/,  this  means  array  errors  of  x  dB  correspond 
to  an  (■  of  -^10  -Approximate  gain  and  phase  error  equivalents  of  various  arra,\'-error  powers 

are  given  in  Table  2- 1 . ' 


2.1.3  Numerical  Implementation  of  the  Bound  Calculations 

W’e  computefl  the  integrated  bounds  given  by  Equation  (1.67)  with  the  derived  in  Sections 
1.3.2  and  1.3.3.  For  evaluating  the  generic-signal  bound,  two  classes  of  distribution  functions  on 
the  signals  were  considered;  we  shall  term  the.se  the  random  constant-envelope  signal  model  and 
the  Gaussian-signal  model.  In  the  random  const  ant -envelope  signal  model,  the  signal  wa\eforms 
were  modeled  as  having  phases  tiiat  were  independent  (iook-to-iook  and  signaJ-to-signal)  random 
variables  uniformly  distributed  over  the  interval  jO. 27r)  with  a  constant  magnitude  for  each  signal. 

’  riie  values  in  this  taljle  w<'r(>  obtained  from  a  Monte  Carlo  simulation  using  10'  trials  for  each 
error  [)ow('r. 
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TABLE  2-1. 


Approximate  Standard  Deviations  of  Gain  and  Phase  Errors  Induced  in  Each  Antenna 
Element  by  Complex  Gaussian  Errors  of  Various  Powers 


Array-Error 
Power (dB) 

Std  Dev.  of 
Gain  Error  (dB) 

Std.  Dev.  of 
Phase  Error  (  ) 

-30 

1.9  X  10  ' 

1.3  i 

-25 

3. .5  V  10' 

2.3 

-20 

0.2  -  10' 

4  1 

-15 

' 

1.1  s  UC 

7.3 

In  the  (laussian-signal  model,  tlie  signal  waveforms  were  modeled  as  being  independent  (look-to- 
look  and  signal-to-signal)  complex  Claussian  random  variables  witli.  possibly  different,  fixed  powers 
for  each  signal.  For  evaluating  the  const  ant -envelujte  signal  bound,  we  considered  only  the  random 
constant-envelope  signal  model. 

The  Fisher  information  matrix  for  the  generic-signal  case  was  computed  using  Equations 
(A. 22).  (.A..'f(J).  (.4. .31).  (A. 38).  (.4.39).  and  (.4.40).  The  matrix  .\y'  was  evaluated  using  Equa¬ 
tion  {1.3G)  or  Equation  (1.37).  as  appropriate.  Note  that  evaluation  of  .\ y'  using  Equation  (1.3G) 
requires  the  inversion  of  only  an  S  ■  5  matrix,  not  an  .V  .V  matrix.  The  information  matrix  for 
the  constant-envelope  signal  case  was  computed  from  the  information  matrix  for  the  generic-signal 
ca.se  using  Equations  (B..3).  (B.7).  and  (B.9). 

The  C'ramer-Rao  bounds  on  DF  for  use  in  Eejnation  (1.G7)  are  given  by  the  diagonal  elements  of 
the  upper-left  5  -  S  submatrix  of  F“'  (see  Section  1.3.1).  This  submatrix  was  computea  according 
to  the  procedur('  outlined  in  .Appendix 

The  expectatifui  integral  in  the  expression  for  the  integrated  bound  iEquation  (1.G7)  was 
approximatefl  by  a  .Monte  Carlo  method  using  3.  100.  1000.  and  1000  Monte  Carlo  trials  for 
experiments  with  100.  10.  2.  and  1  look(s).  respectively.  For  each  Monte  Carlo  trial,  a  sample 
value  oi  the  appropriate  random  siguttl  was  obtained  using  a  pseudorandom  number  generator, 
ami  the  Crtititer-Bao  bound  for  that  .sample  value  was  computed  using  the  procedure  described 
.above.  The  integrated  bound  was  approximated  by  the  mean  of  the  Cramer-Rao  bounds  for  tln' 
signals  genertued.  ( 'fhe  imniber  of  .Monte  Carlo  trials  was  varied  as  a  function  of  looks  since  the 
bound  took  longer  to  compute  when  the  number  of  looks  became  large  and  it  was  found  tiiat  as  the 
number  of  looks  b-  catne  large,  the  bound  became  relatively  independent  of  the  particular  signtd 
waveforms.)  \\ e  al  o  computed  the  standard  error  of  the  estimate  (SEE)  of  this  mean  (i.e..  the  ratio 


of  the  standard  deviation  of  the  Crainer-Rao  bounds  for  the  signals  generated  to  the  square  root 
of  the  number  of  Monte  Carlo  trials)  as  an  index  of  the  error  in  the  approximation.  We  expressed 
the  SEE  as  a  percentage  of  the  mean. 


2.2  RESULTS 

For  scenarios  including  thermal  noise,  we  evaluated  bounds  for  the  following  scenario  parame¬ 
ters:  two  signals,  ID  and  100  looks,  signal-of-interest  (SOI)  ASNR's  of  10.  30.  and  .50  dB  relative  to 
thermal  noise,  iiitcrferer  .4S.\R's  of  10.  30.  and  50  dB.  direction-dependent  array  error  powers  of 
-oc.  -30.  -25.  -20.  and  -15  dB  and  signal  separations  of  0.01.  0.0316.  0.1.  0.316.  and  1.0  B\\  . 

We  also  computed  bounds  on  DF  with  direction-dependent  errors  in  the  ab..Ciice  of  thermal 
noise.  We  shall  refer  to  these  bounds  as  array-errors-only  bounds.  Since,  m  the  absence  of  thermal 
noise,  the  probability  distribution  function  of  the  observation  is  singular  for  .V  >  5  (where  A’  is 
the  number  of  looks  and  S  is  the  number  of  sienaL^.  we  restricted  ourselves  to  scenarios  where 
A'  =  5  when  thermal  noise  was  absent.  (C'ramer-Rao  bounds  can  be  generalized  to  handle  singular 
distributions,  such  as  those  that  oc.jur  in  the  N  >  S  case  here.  Moreover,  when  this  is  done  for 
generic  signals,  the  array-eriors-only  bound  in  the  N  >  S  case  can  be  shown  to  be  equal  to  the 
analogous  bound  in  the  A’  =  5  case.  However,  the  necessary  analysis  would  take  us  too  far  afield, 
and  so  is  no*^  included  here.) 

i\Ionte  Carlo  estimation  of  the  expectation  of  the  Cramer-Rao  bound  (see  Section  2.1.3)  was 
found  to  be  quite  satisfactory,  with  an  SEE  of  the  mean  of  the  bound  on  variance  of  less  than  9% 
in  all  ca.ses.  This  corresponds  approximately  to  an  SEE  of  the  mean  of  the  bound  on  standard 
deviation  of  less  than  5%.  (To  see  this,  note  that  a  small  fractional  change  a  in  x  results  in 
approximately  a  fractional  change  of  2a  in  x^.  since 


(t  ±  ax)^  =  x^±2ax^  +  a‘x^ 

~  x“  ±  2ax~.)  (2.11 ) 

When  the  ASNR  and/or  number  of  looks  was  high,  the  SEE  of  the  mean  of  the  bound  was  much 
smaller. 

2.2.1  Array-Errors-Only  Bounds 

In  the  one-signal,  one-look  case,  it  is  easily  shown  from  Equations  (1.40)  and  (1.48)  that  the 
models,  and  hence  the  Cramer-Rao  bounds  for  the  unknown  generic  and  constant -envelope  signals, 
are  identical.  This  also  holds  in  the  absence  of  thermal  noise,  when  the  only  DF  perturbations 
are  caused  by  the  array  errors.  In  this  case,  the  observation  simply  consists  of  an  unknown  scalar 
multiple  of  the  perturbed  array-response  vector.  It  is  shown  in  Appendix  D  that  this  bound  is 
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LOWER  BOUNDS  ON  DF 


•1  CE  SIGNAL  BOUND 
RANDOM  CE  SIGNALS 
UNIFORM  LINEAR  ARRAY 
10  ANTENNA  ELEMENTS 
1  SIGNAL 
1  LOOK 

-  NO  THERMAL  NOISE 

-  -30  dB  ARRA  Y  ERRORS 

-25  dB  ARRAY  ERRORS 

-J  -20  dB  ARRA  Y  ERRORS 

10° _ 

-15  dB  ARRAY  ERRORS 

Figure  2-2.  Generic  or  constant-envelope  signal  Cramer-Rao  bounds  -  1  signal.  1  look, 
no  thermal  noise. 

independent  of  the  signal  waveform  sample.  The  results  of  numerically  evaluating  this  Cramer-Rao 
bound  are  presented  in  Figure  2-2. 

In  the  two-signal,  two-look  case,  the  Cramer-Rao  bound  on  DF  of  generic  signals  with  array 
errors  only  was  found  numerically  to  be  independent  of  the  sample  values  of  the  signal  waveforms, 
depending  only  on  the  power  of  the  array  errors  and  the  signal  separation.  This  bound  is  plot¬ 
ted  in  Figure  2-3.  The  bound  increased  with  the  power  of  the  array  errors,  was  largest  at  the 
0.316  BW  signal  separation,  and  approached  the  one-signal,  array-errors-only  bound  as  the  signal 
separation  approached  0  (c.f..  Figure  2-2).  In  contrast,  the  array-errois-only  Cramer-Rao  bound 
for  constant-envelope  signals  was  found  to  depend  on  the  signal  waveforms.  The  results  of  nu¬ 
merically  evaluating  the  integrated  version  of  this  bound  with  random  constant-envelope  signals 
arc  presented  in  Figure  2-4.  The  bound  attained  a  maximum  at  moderate  signal  separations  (0.1 
and  0.316  BW)  and  approached  the  one-signal,  array-errors-only  Cramer-Rao  bound  as  the  signal 
separation  approached  0. 
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LOWER  BOUNDS  ON  DF 


GENERIC  SIGNAL  BOUND 
RANDOM  CE  SIGNALS 
UNIFORM  LINEAR  ARRAY 
10  ANTENNA  ELEMENTS 
2  SIGNALS 
2  LOOKS 

NO  THERMAL  NOISE 

-30  dB  ARRA  Y  ERRORS 
-25  dB  ARRAY  ERRORS 
-20  dB  ARRA  Y  ERRORS 

10  ^  10  '  10°  -1 5  dB  ARRAY  ERRORS 

SOURCE  SEPARATION  (BW) 


Figure  2-3.  Generic  signal  Cramer-Rao  bounds  -  2  signals,  2  looks,  no  thermal  noise. 


2.2.2  Integrated  Generic-Signal  Bound  with  Thermal  Noise 
2.2.2. 1  Evaluation  Using  Random  Constant-Envelope  Signals 

The  results  for  this  case  are  presented  in  Figures  E-1  to  E-18.  They  may  be  summarized  as 
follows: 

Effect  of  signal  separation  In  the  absence  of  array  errors,  the  bound  increased  monoton- 
ically  and  unboundedly  as  the  signal  separation  decreased.  When  array  errors  were  included,  the 
bound  for  signals  at  large  separations  generally  increased,  particularly  for  high  SOI  ASNR  and  a 
large  number  of  looks.  However,  tlie  bound  for  signals  at  small  separations  generally  decreased.  In 
addition,  for  a  fixed  power  of  array  error,  the  bound  became  wel'  behaved  as  the  signal  separation 
decrea;sed.  With  all  of  the  array-error  powers  tested,  the  bound  generally  was  largest  at  moderate 
signal  separations  and  smallest  at  the  extreme  separations  of  0.01  and  1.0  BW. 

Effect  of  array-error  power  With  high  SOI  ASNR  and  numbers  of  looks,  the  bound 
decreased  as  the  (nonzero)  power  of  the  array  errors  decreased.  However,  when  the  SOI  ASNR  and 
the  number  of  looks  were  sufficientlj’  small,  the  opposite  occurred. 
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LOWER  BOUNDS  ON  DF 


CE  SIGNAL  BOUND 
RANDOM  CE  SIGNALS 
UNIFORM  LINEAR  ARRAY 
10  ANTENNA  ELEMENTS 
2  SIGNALS 
2  LOOKS 

NO  THERMAL  NOISE 

-30  dB  ARRA  Y  ERRORS 
-25  dB  ARRA  Y  ERRORS 
-20  dB  ARRA  Y  ERRORS 
-15  dB  'array  'errors 

SOURCE  SEPARATION  (BW) 

Figure  2-4.  Constant-envelope  signal  bounds  -  2  random  constant-envelope  signals.  2 
looks,  no  thermal  noise. 


Effects  of  signal  power  and  number  of  looks  In  the  absence  of  array  errors,  the  bound 
decreased  as  the  power  of  the  SOI  or  the  number  of  looks  increased.  It  wcis  independent  of  the  power 
of  the  interferer.  In  the  presence  of  array  errors,  the  bound  continued  to  decrease  as  the  power  of 
the  SOI  or  the  number  of  looks  increased  when  the  SOI  ASNR  and  number  of  looks  were  small. 
However,  given  a  sufficient  SOI  ASNR  and  number  of  looks,  the  bound  became  almost  independent 
of  number  of  looks  and  signal  powers,  numerically  approaching  the  two-signal,  array-errors-only 
Cramer-Rao  bound.  The  bound  for  a  signal  at  a  given  power  tended  to  decrease  to  a  limiting  value 
when  the  power  of  the  other  signal  was  increased.  The  tendency  to  decrease  was  most  noticeable 
when  the  signal  separations,  the  ASNR  of  both  signals,  the  number  of  looks,  and  the  array-error 
powers  were  small. 


2. 2. 2. 2  Evaluation  Using  Gaussian  Signals 

The  results  obtained  by  evaluating  the  generic-signal  bound  with  Gaussian  rather  than  random 
constant -envelope  signals  are  presented  in  Figures  E-19  to  Ei-36.  These  results,  though  slighter 
higher  in  the  10-look  ca,ses.  are  qualitatively  similar  to  the  results  for  random  constant-envelope 
signals. 
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2.2.3  Integrated  Unknown  Constant-Envelope  Signal  Bound  with  Thermal  Noise 

The  results  for  this  case  are  presented  in  Figures  El-37  to  E-54.  They  may  be  summarized  as 
follows, 

2.2.3. 1  Effect  of  Signal  Separation 

In  the  absence  of  array  errors,  the  bound  ha^  a  very  weak  dependence  on  signal  separation 
with  a  maximum  at  moderate  signal  separations.  This  pattern  persisted  in  the  presence  of  array 
errors,  with  the  boiipH  V'ceoming  numerically  independent  of  signal  separation,  given  sufficient  SOI 
ASNR  and  number  of  looks.  Both  with  and  without  array  errors,  the  constant-envelope  signal 
hounds  were  dose  to  the  analogous  generic-signal  bounds  at  a  sigital  separation  of  one  B\V. 

2. 2. 3. 2  Effect  of  Array-Error  Power 

Unlike  the  generic-signal  bound,  the  constant-envelope  signal  bound  generally  increa.sed  with 
increasing  array-error  power  at  all  SOI  ASNR's  and  numbers  of  looks. 

2. 2. 3. 3  Effects  of  Signal  Power  and  Number  of  Looks 

In  the  absence  of  array  errors,  the  bound  decreased  with  increasing  power  of  the  SOI  or  number 
of  looks.  In  the  presence  of  array  errors,  a  similar  pattern  occurred  when  the  SOI  ASNR  and 
number  of  looks  were  small.  However,  given  sufficient  SOI  ASNR  and  number  of  looks,  the  bound 
numerically  approached  the  one-signal,  array-errors-only  Cramer-Rao  bound.  Thus,  in  contrast  to 
the  limiting  bound  for  the  generic-signal  case,  the  limiting  bound  for  the  constant -envelope  case 
was  independent  of  signal  separation.  The  bounds  (with  and  without  array  errors)  showed  little 
variation  with  the  power  of  the  interferer.  The  constant-envelope  signal  bounds  tended  to  be  more 
optimistic  than  the  generic-signal  bounds.  When  array  errors  were  present,  this  difference  was  most 
noticeable  at  moderate  signal  separations,  when  the  SOI  ASNR  and  number  of  looks  were  low  and 
the  array  errors  were  small.  When  array  errors  were  absent,  this  difference  was  largest  at  small 
signal  separations  and  quite  significant  at  all  ASNRs  and  numbers  of  looks  examined. 

2.3  DISCUSSION 

Several  aspects  of  the  results  bear  some  discussion.  One  significant  aspect  is  the  tendency  of 
the  bound  on  DF  of  a  generic  signal  to  decrease  as  the  power  of  the  generic  interferer  increases. 
1  his  may  be  due  to  the  fact  that  the  subspace  spanned  by  the  array- response  vectors  (i.e..  the  signal 
subspace)  becomes  easier  to  estimate  as  the  power  of  the  interferer  increases.  While  the  bound 
is  not  necessarily  achievable,  it  suggests  that  the  presence  of  direction-dependent  errors  docs  not 
preclude  accurate  DF  of  signals  at  low  signal-to-interference  power  ratios. 


Aiiotlier  significant  aspect  is  that  the  presence  of  small  array  errors  improved  the  bound  on  DF 
of  generic  signals  at  small  signal  separations.  This  finding  may  be  due,  in  part,  to  the  fact  that  array 
errors  increase  the  mean  distance  (in  complex  A/-space)  between  the  actual  array-response  vectors 
of  the  signals  and  thus  result  in  an  apparent  increase  in  signal  separation.  However,  the  increa.se 
in  mean  distance  between  the  array-response  vectors  does  not  entirely  explain  this  phenomenon, 
since  the  bound  with  arra>-  errors  actually  decreased  sigtdficantly  as  the  signal  separation  went 
from  moderate  to  small.  While  the  bound  is  again  not  necessarily’  achievable,  it  is  consistent  with 
the  possibility  of  performing  DF  with  direction-dependent  enurs  and  small  signal  separations. 

In  the  geneiic-.sigiuil  r-ise.  the  integrated  bound  on  DF  of  2  signals  with  array  errors  and 
thermal  noi.se  approached  the  2-si^nal.  array-errors-only  Cramer-Rao  bound  as  the  SOI  .AS.MR  and 
the  number  of  looks  became  sufficiently  high.  Tiie  array-errors-only  Crame'r-Rao  bound  was  found 
to  be  independent  of  the  signal  waveforms,  depending  only  on  the  power  of  the  array  errors  and 
the  signal  separation.  Tfiis  finding  is  consistent  with  the  notion  that.  gi\'en  a  sufficient  SOI  .\SN'R 
and  number  of  looks,  the  array  errors  are  the  dominant  factor  limiting  DF.  The  array-errors-only 
C'r;uner-f{ao  bound  was  found  tir  be  markedly  different  from  the  integrated  bound  on  DF  with 
tiiermal  mrise.  but  no  array  errors,  with  respect  to  its  dependence  on  signal  separation.  The  latter 
liound  increased  sharply  as  the  signal  separation  decreased,  while  the  former  did  not.  In  fact,  the 
array-errors-only  Cramer-Rao  bound  for  two  signals  approached  the  array-errors-only  Cramer-Rao 
bound  for  a  single  signal  as  the  signal  separation  approached  zero.  This  latter  bound  is  independent 
of  the  signal  waveforms,  depending  only  on  the  power  of  the  array  errors,  and  is  identical  for  generic 
and  constant-envelope  signals  (see  Section  2.2.1). 

With  respect  to  the  integrated  bound  for  constant-envelope  signals,  the  results  are  significant 
in  that,  in  many  ca.ses.  in  spite  of  the  presence  of  array  errors,  the  bound  remained  more  opti¬ 
mistic  than  the  analogous  bound  for  generic  signals,  particularly  in  scenarios  with  moderate  signal 
separations,  small  array  errors,  and  a  low  SOI  ASNR  and  number  of  looks.  This  suggests  some 
potential  improvement  in  DF  still  exists  through  use  of  the  constant-envelope  signal  assumption 
in  the  presence  of  small  direction-dependent  errors,  though  not  to  the  extent  that  is  true  in  the 
absence  of  errors. 

In  the  constant-envelope  signal  case,  the  integrated  bound  on  DF  of  2  signals  with  array  errors 
and  thermal  noise  approached  the  single-signal,  array-errors-only  Cramer-Rao  bound  as  the  SOI 
.\.S.\R  and  number  of  looks  became  sufficiently  high.  The  latter  bound  is  simply  the  bound  on  DF 
of  a  signal  given  an  observation  consisting  of  an  unknown  scalar  multiple  of  the  perturbed  array- 
response  vector  of  that  signal  (see  Section  2.2.1).  This  finding  is  consistent  with  the  notion  that, 
given  a  sufficient  SOI  .ASNR  and  number  of  looks,  the  array-response  vector  of  the  SOI  becomes 
determinable  to  within  a  complex  scalar.  Since  the  dependence  of  the  olrservation  on  the  direction 
occurs  ordy  through  this  vector  and  the  array  errors  on  this  vector  were  modeled  as  independent 
of  all  the  other  random  quantities  in  tlie  problem,  we  may  conclude  that  the  interfering  signal  lias 
little  degrading  imjtact  on  DF. 
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3.  CONCLUSIONS 


C  rarner-Rao  lower  bounds,  on  the  variance  of  unbiased  DF  algorithms,  were  derived  for  un¬ 
known  generic  and  constant-envelope  signals  received  with  complex  Gaussian  direction-dependent 
array  errors  and  thermal  noise.  These  bounds  were,  in  general,  functions  of  the  specific  complex 
amplitudes  of  the  signals.  The  Cramer-Rao  bounds  w'ere  used  to  obtain  other  lower  bounds,  which 
we  term  integrated  bounds,  on  the  variance  of  DF  algorithms  for  th*^  ease  where  the  complex  am¬ 
plitudes  of  the  signals  have  a  probabilistic  distribution  which  is  unknown  to  the  algorithm.  The 
integrated  bounds  have  the  desirable  property  of  not  depending  on  the  specific  complex  amplitudes 
of  the  signals.  In  some  cases,  the  Cramer-Rao  bounds  on  DF  turn  out  to  be  functionally  independent 
of  the  specific  complex  amplitudes  of  the  signals.  In  such  cases,  the  integrated  bound  is  numerically- 
identical  to  the  Crarner-Rao  bound  which  is  obtained  for  any  allowable  set  of  complex  amplitudes. 
We  evaluated  the  bounds  under  various  conditions,  for  up  to  two  generic  or  constant-envelope  sig¬ 
nals.  Those  Cramer-Rao  bounds  which  were  independent  of  the  specific  complex  amplitudes  of  the 
signals  were  evaluated  exactly,  while  the  integrated  bounds  were  evaluated  using  a  Monte  Carlo 
simulat.jn  to  approximate  expectation  integrals  with  respect  to  the  probability  distributions  of  the 
complex  amplitudes  [see  Equation  (1.67)j.  The  behavior  of  the  bounds  is  briefly  summarized  below; 
notizero  array-error  powers  are  assumed  unless  otherwise  noted. 

1.  I'nknown  generic-signal  bound: 

(a)  The  integrated  bound  tends  to  decrea.se  with  increasing 
interferer  power,  approaching  a  limiting  value  as  the  power 
of  the  interferer  gets  large. 

(b)  At  small  signal  separations,  the  integrated  bound  without 
array  errors  increases  unboundedly  as  the  signal  separation 
decreases.  The  addition  of  small  array  errors  results  in  a 
marked  decrease  in  the  integrated  bound  at  small  signal 
separations,  with  the  bound  becoming  well  behaved  as  the 
signal  separation  decreases.  These  results  suggest  that  the 
presence  of  array  errors  may  actually  improve  DF  of  closely- 
spaced  signals. 

(cj  With  a  sufficiently  large  SOI  ASNR  and  number  of 
looks,  the  integrated  bound  for  two  signals  approaches  the 
Cramer-Rao  bound  for  DF  of  two  signals  (given  two  looks) 
with  array  errors  only.  The  latter  bound  was  numerically 
independent  of  the  complex  amplitudes  of  the  signals. 

(d)  The  array-errors-only  bound  on  DF  of  two  signals  (given 
two  looks)  depends  only  on  the  power  of  the  array  errors 
and  the  signal  .separation.  The  bound  increa.ses  with  the 
power  of  the  array  errors,  is  largest  at  moderate  signal 
separations,  and  approaches  the  single-signal  (one-look) 


bound  as  the  signal  separation  approaches  0.  The  latter 
bound  is  also  independent  of  the  signal  waveform. 


2.  Unknown  constant-envelope  signal  bound: 

(a)  The  integrated  bounds,  both  with  and  without  array  er¬ 
rors,  are  only  weakly  dependent  on  the  power  of  the  inter- 
ferer  and  signal  separation,  and  are  close  to  the  analogous 
bounds  for  generic  signals  at  a  signal  separation  of  one  B\V. 

(b)  With  no  array  errors,  the  integrated  bound  for  constant- 
envelope  signals  is  far  more  optimistic  at  small  signal  sep¬ 
arations  than  the  integrated  bound  for  generic  signals,  as 
can  be  concluded  from  observations  (lb)  and  (2a)  above. 
With  array  errors,  the  integrated  bound  for  constant- 
envelope  signals  is  still  somewhat  more  optimistic  than  the 
analogous  bound  for  generic  signals,  particularly  at  small 
array  errors,  low  SOI  ASNR  and  numbers  of  looks,  and 
moderate  signal  separations. 

(c)  With  a  sufficiently  large  SOI  ASNR  and  number  of 
looks,  the  integrated  bound  for  two  signals  approaches  the 
Cramer-Rao  bound  for  DF  of  one  signal  with  array  errors 
only.  The  latter  bound  is  identical  to  that  for  a  single 
generic  signal,  discussed  in  (Id)  above. 
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4.  FUTURE  WORK 


Ft)r  practical  applicatioiiN,  it  would  be  useful  to  liave  bounds  vcith  more  realistic  models  for  the 
direc’tioti-depemient  errors.  For  example,  one  should  allow  for  errors  of  unknown,  possibly  difierent 
variances  fur  each  sigmtl  with  irossible  element-to-element  and  or  signal-to-signal  error  correlations. 
One  should  also  incorporate  other  error  sources  into  the  model,  such  as  multiplicative  direction- 
independent  errors  and  de\  iatiotis  of  th<‘  signals  from  perfect  constant-envelopeness.  Finally,  arrays 
designed  to  have  dillerent  response  vectors  for.  say,  vertically  vs.  horizontally  polarized  sources  could 
hiive  ditferent  sensit i\  it ies  to  the  various  error  sources,  in  terms  of  their  DF  capabilities. 


APPENDIX  A 

THE  INFORMATION  MATRIX  FOR  GENERIC  SIGNALS 


In  the  derivation  that  follows,  we  shall  use  a  number  of  defined  terms.  For  the  convenience  of 
the  reader,  we  list  a  number  of  these  here. 

We  define  the  F  >  F  matrices  IF,  li'.  and  IF  by 

IF^^-^F,  (A.F) 


\V'^V><V. 


1'F=^'-^^F  (A.. FI 

(the  noltitions  U  and  ll  are  mnemonic  only;  they  do  not  indicate  derivatives  of  IF),  the  S  ■  .S' 
Hennitian  matrix  /’  by 

P  (A,4) 


the  S  .'V  matrix  Q  by 


the  5- vector  by 


-  dyif  r'l 

Qn  ^-Ti¬ 


the  .S'  .S'  matrices  D„  by 

2{Cr,ql)^\V\ 

the  S  y  .S’  Hennitian  matrix  B  bv 


B  IF  +  M(^P\ 


the  .S'  '  .S'  matrices  by 

2(A  v' )„„'(<••; 

and  the  .S’  ■  .S'  matrices  Hm,'  by 


(A. 10) 
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When  (T'  >  0.  F  is  given  explicitly  by 

P  =  -  ^C^Cr'C^'iC  i  En,  1.36 


(T-  cr- 

\Mien  fT-  =1  0.  A’  =  5.  and  (~  >  0.  P  is  given  explicitly  by 

P  =  j  3-  ) 


-Es- 


We  will  now  dcri\'e  an  explicit  expression  for  the  Fisher  inlorniation  matrix  in  the 
case.  The  derivatives  in  Equation  (1-62)  are  given  by 

ar,  or,  or, 


and 


^  _  2/^^V// 

(tr,  dr, 


dC^' 

1m 

or. 


=  ( 


H 


V.  I 


M- 


Evaluating  these  equations  for  the  various  components  of  f,  one  obtains 
dfi 


d'lf 


=  (C  ■  V'£.„)l5  (12.  Eq.  E291). 


dpi 

C)(l„, 


( 


dc 


ddfis 


y)is 


dCnfi 

dans 


(Ens 


V)h- 


djl^ 

do  ns 


,  dC 
dOns 


riFs- 


'p^iEns  ■  V  IFs. 

dOns 


(A, 11) 

(A.  12) 
generic-signal 

(A.t3) 

(A.14) 

(A.15) 

(A.16) 

(A.  17) 
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(A.IM 


()a„s  Vwi„.s 


'  ""■''  \da„J 


CE,,  ■:/.„. 


(A. 19) 


=  r\p^i:.,.c‘ 


W'c  shall  now  evahiatc  the  F starting  with  Fu,,-.  Substituting  into  Equation  (l.()2) 


^  «  {|£ 


-  2Efr{/{F  •  V-E,,)^'(Av'  >:  I.mHC  • 


2R  l‘J'(C^' .\;^-^C)Z{EilV^V E,-s'y\s\  (Eq.  1,18) 


2R:r" 


P- (  F"  11  ■  F,,ee ),  1  ,s  1  ( Eqs.  A  .4  and  A .3 ) 


-  21? ir(/:pz(ir).yF,,qrs|  (:2.  Eq.  1.32;) 

-  2h'  [i'JtP),,'{U')ss'E,.,Ps] 

^  2i?[(F),,dir)«'j . 


21?(/’Air), 


Next  we  address  tiu’  /•  •  ..  wiierc'  eriiials  (7„  or  o„.  Define  b„s  and  h„  by 


,  d.  f  ?^e„, 

F,.v  - 


(A.23) 
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Thus  for  .i,i  = 

h„  =  L-,i  (Eqs.  1.40  and  1.45). 
wliile  for  .1,,  =  o„. 

b„  =  ja„^c„  (Eqs.  1.40  and  1.45). 


Substituting  into  Equation  (1.62)  gives 


(F 


=  2?K 


/  dp 


I  i  JliL 

V  c)3ns  ) 


=  2']i[b:JF(E,„  ■  r)"(Av‘  ;■  /,/) 

x(C'  .  r£'.,<s/)r5j  (Eqs.  A. 15.  A. 17  and  A. 16) 

=  2Ji\b:J^i(E,„\:^^C)Z}(V^VEp,-)]ls]  (Eq.  1.18) 

=  2'}?  ^:J^l(EsnQ^)^{\VEs's')3s]  (Eqs.  A.2  and  A. 5) 
=  2'^[b’J^(Q^)nA^V)ss'Ess-ls]  (Eq.  1.13) 

=  21?  [6;„((?")„,-(ir),,-] 

=  2^[b:jQ)A{^Vhs] 

=  2=R{V«(Q)..'n:*(»E)..,'} 

=  21?  {  I.-; ‘Jli ■)...«'}  (Eq.  A.6) 

2 1?{f6;r///),,dH  •).,.«-} 

=  2’R{{b„ql,)...A'^y‘)ss'}. 

Thus 

-  2f?:(6„qt)-lE‘;. 

Substituting  for  the  3„.  one  obtains 

F^„,  -  AR'xCn<g)GWA  (Eq.  A.25) 


=  mOr,)  (Eq.  A.7) 


(A. 24) 


(A.25) 


(A.26) 


(A.27) 


(A.2S) 

(A.  29) 


(A. 30) 
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1 


aiul 


F 


Lastly  we  address  the  Fj 
Define  and  gn  by 

cief 

”  d- 

!n,< 


2;R{;Oa„F<.'„)9l'pU-*}  (Eq.  A.26) 

-2:s{  {an^V„)qlp\'V'} 

-2A{:(d„=P)(ii:c^  )!□»•'} 

-2A{;(d„nP)(r5nq;,ppii-‘}  (Eq.  1. 11) 

-2.p(PlP)-:(P(fp3li-*}  (Eq.  1.12) 

-(d„p)Lj2A[(Pq^pU-j  (Eq.  1.10) 

(Eq.  A.7).  (A.31 

where  is  equal  to  a,,  or  o„  and  P'  is  equal  to  d^>  or  o,,' 


and 


-  def  r 

S>i  ~  ffn  I  9n2 


ffuS  j  • 


(A.32) 


(A.33) 


Substituting  into  Equation  (1.62)  gives 


(f =  21? 

■Jri 

Tlie  left-hand  term  evaluates  to 

27m 


Oft 

03,,.. 


dg 


y 


dg 


d 


in'  s’ 


^tr(  A-’:^A-'  . 


(A.34) 


03,, 

J 

2:« 

^^n.s.Fi'.e 

r,"(E„ 

21?  Ih*  r/.s-. 

^  ua‘2'i 

dslF, 

2«  ! 

i 

21? 

[iKSl) 

-'(Ap 

(21?  {(Ap 

3F 


-I  . 


;(F,,^A7.‘F„-,0a(r"V)jrs]  (Eq.  1.18) 


(A.  35) 
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The  right-liand  term  evaluates  to 


t)\  dA 


tr 


=  tr  {(A  v'  X /a/)>‘ 


o  r  dCri^f' 


:  2  [ 


F 


X  r  ^-^En's'C^  +  i,  /Mi}  (Eq.  A.lh  aud  A.20) 

‘  Lt/'n's'  VcE-n's'/ 

=  (Ux{[XlHhnsE,,sC^ 

y \^^^[gn's'E„>_,'C^  -r  g*„,^,CE^,„')]  s  hi}  (Eq.  1.7) 

=  tr  i(6,„5„-,-A^.'£„,(:’"A;^.’E„,,-C"  ^  6„,5;,,,A  v’£.„C"A  v'CE,-,,. 

-^l,^/r^^eAv'^£.,MVv  +  6;.5:-.-Av'r£,„A^-'C£,.„0  •;  I.m] 

=  M(^  tr  (6,,.9,.-,-A;;.'£,,r"A;^.'£„-,-C"  +  bn,g'„,,,A:^^E„sC^\:^^CE,>r,' 

^bh9n's'A:xCEs„.\:^^E„>,,C^  +  E„X's'-^:xCE,r,.\s'CE,.„')  ( '2.  p.30;) 

=  Me^  tr  (b.sgn’s'Kx  Er,,C^ E„>s'C^  -  6,„5;'.'A,v'£nsC"A;'C£,-„- 
-bhgn's'En’s'C^A:^^CEsnA:^:^  +  6;,5:'.-C£yn'A)^.'C£,.A^.’ )  (Eq.  1.1 ) 

=  2.\h^'^[tr  [bns9n's'A-^^E„sC^Ay^En','C^ 

■^b:,,g,>,'En's'C^A-^^CEsnAi:^)] 

=  2.1/rA'E[tr  (6,„p„-,.£,,C’"A:,-'£n'.xC" A)^.' 

-^^;,,9n'.'£„vr"A:,.'r£*nA,v*)]  (Eq.  1.1) 

=  2M(^^'  U  (b„sgn's'EnsQE„'s'Q 

^fj;,gn',,/£„'s'P£,nA  v')j  (Eqs.  A. 4  and  A. 5) 

=  2.Vf'‘5?  ’  tr  {bnsgn'.s'{Q)sn'E„s'Q 

-b:„g„,.{P)s'sE„>„A-^')]  ([2.  Eq.  1.3lj) 

=  2.yf‘‘3?  [6„,,g„.,.<(Q).,„-(Q)s'„  4  6;,,9„'5'(P)s-s(A^vVlnn']  (Eqs.  1.2  and  1.1) 


CE,„  x: /mI(A),.’ :•  /m) 


2Mhji  [(br,gJhs.Aqn'Qn  ).■,.■<'  +  (6;5n'W(P*)ss'(A,vMnn'] 
{2Mh^{[{bng'!,^)a(qr.'gn)]  +  (A^^'hn'iKgl’MPh}) ■  (A.36) 
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Thus  the  -  are  given  bv 

=  23f?{(A;^.iW(6;5lOcir 

+  Me^[{bngj;.)D(q„.ql)  +  )nn'iK9l')°n] 

=  23f?{(Avi  W(6;5X')D(»'  + 

+  A[e'^{bn9l')0{q„'ql)} 

=  25R{(A^.M„„-(6;5r0DS  +  ^^/f'’(^n5?)a(<rn'9t)}  (Eq.  A.8). 

Substituting  for  the  J„  and  one  obtains 

=  2^{{A:^^)„n'{t-:iZ)^B  +  Me\vniZM 

=  5R{Gw  +  ^nn')  (Eqs.  A.9  and  A.IO), 


(A,37) 


(A.3S) 


+Me‘^\{janOil)n)il^>]0{q„,ql )} 

=  25i{-j{(A^.'W[(a„at^n)’^^-lD5 
-A/€"[(a„ar?n)v^<]D{9n'gn)}} 

=  2^{{A-^^)nn'\{a„ntp,,ri^.pB 

-Me*[(a„0i;r,)V^,]a(q„.qj;)} 

=  29{(A;;.'W(a„iI)D(v^;v^-)D5 

-A/£‘’(a;i5)D{v„i^-)0(9n'gD}  (Eq-  1-12) 

=  2(a„ll)a^i{(A;v‘)nn'(t^>T')D5 

-Me\^fn^'Mqn’ql)} 

=  (a„i^)D9(Gw -//„„')  {Eqs.  A.9  and  A.IO),  (A.39) 

^ a  dT  ~  25R{(  Ay^.' )„n'[(ja„Dt/^n) 

+  A/f  ■*  [(  jd„  Dt/Jn  ) ( ja„' Dt^n' )^1 (9n' ?! )  } 

=  2iR{(A;;;)„„-[(a„DV^n)*(an'DtPn.)^l°^ 

-  M('^\ianOp„){a„>Dip„>)'^]0{q„,ql)} 
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-A/e'*(a„a^,)D(v5’„i/^0°(9n'9n)} 

=  2{a„al)n3i{{A],})„r,'{^'r,^')OB 

-Mt‘^{pn'>i^')0{qn'ql)} 

=  (a„o^',)03R(G„„-  -  H„„')  (Eqs.  A.9  and  A. 10). 


(A.40) 
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APPENDIX  B 

THE  INFORMATION  MATRIX  FOR  CONSTANT-ENVELOPE  SIGNALS 


If  one  makes  the  icieutificat itni 


the  expressions  for  Op  0o„<..  d.\  ihif.  and  0.\  Oons  given  by  Equations  (A. 15).  (A. 17). 

( A.lf<).  and  (.A. 20).  respectively,  remain  valid,  since  the  observation  and  its  functional  dependence 
on  u  and  the  o„  remain  unchanged.  W'e  then  only  need  to  evaluate 

Op  OC 

- —  —  •  I  )l.s  (Eq.  .A.i:l) 


=  •  i-iG 


=  E: 


(Etp  A.  10) 


d(y,  \0c\f: ) 


I.M  (Eq.  A.  14) 


=  E- 


=  y,^'p±E.,.c" 


=  X  t—  (Eq.  A.IO).  (B.d) 

where  tiu'  partial  derivatives  with  res[)ect  to  in  Equations  (B.2)  and  (B.3)  are  understood  to 

b('  flelined  as  in  tlu'  generic-signal  case,  to  obtain  expressions  for  all  the  Op  Or,  and  fl.V  Or,. 

Since  the  partial  derivatives  of  p  and  with  respect  to  p  and  the  o„  are  the  same  as  in  the 
generic-signal  case,  it  follows  that  the  with  0  4-  «  and  y  A  o  are  given  by  the  corresiionding 
ex[)ressions  in  .Appendix  .A.  It  remains  to  evaluate 


3  (  ^^0 


ly) 


(Eq.  1.62) 


which  gives 


n  =  1 


A' 


V  Udj^g  / 


-1 

du.1 


(Eq.  B.2) 


Y^{fa„a)s>'  (Eq.  A.27). 


>1  =  1 


(B.4) 


which  gives 


and 


which  gives 


f'as  = 

n  —  \ 


dp  ,-i  dfl 

5q  ,  /  "  da 


tr 


dag  dag' 


x  x 


=  EE2» 

S  =  l7l'=  1 

.V  v 


-1  dp 


da„' 


ST'  f  d-*^  .  - 1  d\ 

(^a„d„,  )s5'- 

r2  =  ln'  =  l 


^aa  -  51  51 

n  =  ln'  =  l 


=  2.« 


(9o 


dp  \"a-^ 


(B.5) 


(Eq.  1.62) 


(B.6) 


(B.7) 


,  aA  ^A 
tr  A  - A  ^ 


d(pns  dag 


(Eq.  1.62) 


.V 


5^23? 


>i'  =  i 
.V 


dp  dp 


da„‘g' 


^  d\  d.\ 

^tr  A  1,; - A 


n'  =  l 


51  ^^o„S„, 

n'  =  1 


.\ 


F-  .  =  y  F-  -  . 


(B.8) 


(B.9) 


n'=  \ 


Thus  the  information  matrix  for  tlie  unknown  constant-envelope  signal  case  is  easily  constructed 
from  the  information  matrix  for  the  generic-signal  case. 
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APPENDIX  C 

NUMERICAL  REDUCTION  OF  THE  INFORMATION  MATRIX 


For  notational  convenience  in  this  section,  let  Um-  denote  the  upper-left  m  x  m  submatrix  of  a 
positive-definite  symmetric  A'  x  A’  matrix  V.  Suppose  F  is  the  L  x  L  Fisher  information  matrix  for 
the  L-dimensional  vector  of  unknown  parameters,  r.  As  derived  in  Section  1.3.1.  the  Cramer- Rao 
bounds  on  DF  for  the  individual  signals  are  given  by  the  diagonal  elements  of  the  upper-left  S  x  S 
submatrix  of  F'*.  (F^^)  5  .  We  shall  now  describe  a  procedure  for  obtaining  (F^^)l5-. 

f’artition  V  as 


Define  the  function  /? v  mapping  llie  set  of  A'  x  A’  positive-definite  symmetric  matrices  to  the  set 
of  (.V  -  1)  X  (A’  —  1)  positive-definite  .symmetric  matrices  by 


:  ( i  )  -  D  A'  - 1 


Thus  we  have  that. 


=  \R^{y)r\  (C.3) 

(Eq.  1.21).  We  shall  show  that  {F“^)iqj  can  be  obtained  by  applying.  L  -  q  times,  the  operator  R 
(appropriately  dimensioned)  to  F,  i.e.,  that 

(F-')[<,;  =  [R,.iR,^2--Rl{F)]-^.  (C.4) 

The  result  is  true  for  g  =  F  -  1  by  Equation  (C.3).  Assume  (the  induction  hypothesis)  that  it 
is  true  for  9  -x  1,  Then 

=  {F,.i!/r:,+2/Z,+3"-Fi(F)]}-'  (Eq.  C.3) 

=  [F,.iV2  •Fl(F)]-’  (C.5) 

(the  first  line  merely  identifies  the  upper-left  q  x  q  submatrix  of  F~^  with  the  upper-left  q  x  q 
submatrix  of  the  upper-left  (9  1)  x  (g  -t-  1)  submatrix  of  F“').  The  result  is  thus  true  for  q.  and 

the  theorem  follows  by  mathematical  induction. 

We  used  the  technique  described  above  to  obtain  (F“')[5i. 
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APPENDIX  D 

THE  ARRAY-ERRORS-ONLY  CRAMER-RAO  BOUND  FOR  5  =  .V  -  1 


In  this  appendix,  we  shall  evaluate  the  array-errors-only  Cramer- Rao  bound  for  the  case  where 
5  =  .V  =  1.  As  noted  in  Section  2.2.1,  the  models  for  generic  and  constant-envelope  signals  are 
identical  for  this  ctise.  Thus,  with  no  loss  of  generality,  we  shall  compute  this  bound  using  the 
notation  of  the  generic-signal  case.  For  simplicity,  we  shall  assume  that  the  ideal  array-response 
vectors  are  obtained  from  an  array  of  omnidirectional  elements  with  pairwise  symmetry  about  some 
phase- reference  point  in  space.  For  such  an  array,  the  ideal  array-response  vectors  are  of  unit  length 
and  have  real  inner  products  with  each  other.  Letting  Ut  denote  the  true  value  of  u.  we  have 

IF  =  (Eq.  A.l) 

=  1.  (D.l) 

IF  =  v^(ut)v{ut)  (Eq.  A.2) 

=  ■^\r^Uui)v{v)]\u=u, 

=  ^{^\v^{Ut)v(u)]\u  =  u.} 

du 

=  (u,)v{ut)} 

I 

=  +  V  (U()F(u,)l 


= 

(  ft)]  lu  =  TI(  1 

liu=U(  1 

=  0. 

(D.2) 

IF  = 

{ut)v{ut) 

(Eq.  A.3) 

-  1:F( 

(D.3) 

•A.v  --  r'afj 

(D.4) 

(E(|,  l.T'li. 


/ 
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(Eq.  D.4), 


(Eq.  A.4). 


Q  = 


(Eq.  A..5), 


(Eq.  A.7). 


(Eq.  A.8), 


f-Cll 


1 


D  =  0 


B  =  1  +  Mr 


G  = 


2  +  2A/f^ 


(‘■a 


11 


(Eq.  A.9). 


H  = 


f‘'a|ie>2a’u 

2M 


(Eq.  A.IO) 


‘11 


F., 


2||u(u()||^ 


(Eq.  A. 22). 


(Eq.  A.. 30). 


F„,,  =  0 


=  0 


(D.6) 


(D.7) 

(D.8) 

(D.9) 

(D.IO) 

(D.ll) 

(D.12) 

(D.13) 

(D.14) 
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(Eq.  A.31). 


2  f  2M€~  2M 

f'aa  -  - 5-0 - 1 - —  (Eq-  A. 38) 

f  a|j  fljj 

2  +  2M('  +  2M(^ 

^5  7y 

f“«n 
2  +  4A/e^ 

O  O  ' 

€“aji 

Foa  =  0 


(Eq.  A. 39), 


and 


E<j)0  ~ 


‘11 


2  +  2A/e“  2M 


■j  2 


U-jj 

2  -r  2iUf2  -  2A/f2 
^2 

f2- 


(Eq.  A.40) 


E  = 


2|ir(u;  )|!' 


0  0 


(D.15) 


(0.16) 


(D.17) 


(D.18) 


(Eq.  1.60).  Thus  the  Cramer- Rao  bound  on  OF  is  given  by 


var  (it)  > 


(D.19) 


which  is  independent  of  the  complex  amplitudes  of  the  signals.  This  is  equal  to  the  thermal-noise- 
only  Cramer- Rao  bound  on  OF  of  a  generic  signal  with  an  ASNR  of  ^  in  the  S  —  N  =  I  case  [1. 
Eq.  2.17:.  where  the  value  of  \\v{ut}\\^  here  is  equal  to  the  quantity  (27rf')^  in  [1],  where  is  the 
mean-squared  distance  of  the  (linear)  array  elements  from  the  phase  center  of  the  array  in  units  of 
wavelengths-squared. 
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APPENDIX  E 
GRAPHICAL  RESULTS 


This  appendix  gives  a  graphical  presentation  of  the  results  discussed  in  Sections  2.2.2  and 
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Figure  E-1.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  10  looks. 
SOI  ASSR  -  10  dB.  interferer  ASNR  —  10  dB. 
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Figure  F-2.  Generic  signal  hounds  -  2  random  constant-envelope  signals.  10  looks. 
SOI  ASNR  -  10  dB.  interferer  ASNR  —  30  dB 
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Figure  E-5.  Generic  signal  bounds  -  2  random  const  ant -envelope  signals.  10  looks. 
SOI  ASNR  =  30  dB.  interferer  ASNR  =  30  dB. 
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Figure  E-6.  Generic  .signal  hounds  -  2  random  constant -envelope  signals.  10  looks. 
SOI  .ASNR  =  30  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-7.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  W  looks. 
SOI  ASSR  —  50  dB.  interferer  ASKR  ~  10  dB. 
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Figure  E-8.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  lU  looks. 
SOI  ASNR  -  50  dB.  interferer  ASNR  =  30  dB. 
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Figure  E-9.  Generic  signal  bounds  -  2  random  constant-enveiope  signals,  JO  looks. 
SOI  ASNR  =  50  dB,  interferer  ASNR  =  50  dB. 
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Figure  E-10.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  100  looks. 
SOI  ASNR  =  10  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-ll.  Generic  signal  bounds  -  2  random  constant -envelope  signals.  100  looks. 
SOI  ASSR  =  10  dB.  interferer  ASNR  =  30  dB. 


LOWER  BOUNDS  ON  DF 


GENERIC  SIGNAL  BOUND 
RANDOM  CE  SIGNALS 
UNIFORM  LINEAR  ARRAY 
10  ANTENNA  ELEMENTS 
2  SIGNALS 
100  LOOKS 
Pi  =  10  dB  ASNR 
P2  =  50  dB  ASNR 

NO  ARRAY  ERRORS 
-30  dB  ARRA  Y  ERRORS 
^5  dB  ARRA Y ERRORS 
-20  dB' ARRAY  ERRORS 
-15  dB  ARRA  Y  ERRORS 


Figure  E-12.  Generic  signal  bounds  -  2  random  constant -envelope  signals.  100  looks. 
SOI  ASSR  -  10  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-13.  Generic  signal  bounds  -  2  random  constant -envelope  signals,  100  looks, 
SOI  ASNR  =  30  dB,  interfcrer  ASNR  =  10  dB. 
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Figure  E-14.  Generic  signal  bounds  -  2  random  constant-envelope  signals,  100  looks, 
SOI  ASNR  =  30  dB,  interferer  ASNR  =  30  dB 
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Figure  E-15.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  100  looks. 
SOI  ASNR  —  30  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-16. 
SOI ASNR  = 


Generic  signal  hounds  -  2  random  const  ant -envelope  signals.  100  looks. 
50  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-17.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  100  looks. 
SOI  ASNR  =  50  dB.  interferer  ASNR  =  30  dB. 
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Figure  E-18.  Generic  signal  bounds  -  2  random  constant-envelope  signals.  100  looks. 
SOI  ASNR  =  50  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-19.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks.  SOI 
AS.\R  =  10  dB,  interferer  ASNR  =  10  dB. 
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Figure  E-20.  Generic  signal  hounds  -  2  complex  Gaussian  signals.  10  looks,  SOI 
ASNR  =  10  dB.  interferer  ASNR  =  30  dB. 
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Figure  E-2L  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks.  SOI 
ASNR  =  10  dB,  interferer  ASNR  =  50  dB. 
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Figure  E-22.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks.  SOI 
ASNR  =  30  dB,  interferer  ASNR  =  10  dB. 
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Figure  E-23.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks.  SOI 
ASSR  =  30  dB.  interferer  ASNR  —  30  dB. 
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Figure  E-24.  Generic  signal  hounds  -  2  complex  Gaussian  signals.  10  looks,  SOI 
ASNR  -  30  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-25.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks.  SOI 
ASSR  =  50  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-2().  Generic  signal  hounds  -  2  complex  Gaus.sian  signals.  10  looks,  SOI 
ASNR  —  50  dB.  interferer  ASNR  =  30  dB. 
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Figure  E~27.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  10  looks,  SOI 
ASNR  —  50  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-28.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  100  looks.  SOI 
ASNR  =  10  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-33.  (generic  signal  hounds  -  2  complex  Gaussian  signals.  100  looks.  SOI 
AS.\H  =  30  dB,  interferer  .ASNR  =  50  dB. 
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Figure  F,-3I.  Generic  signal  bounds  -  2  complex  Gau.s.sian  signals.  100  1  !:s.  SOI 

ASNR  -  50  dB.  interferer  ASNR  -  10  dli 
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Figure  E-35.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  100  looks.  SOI 
ASNR  =  50  dB.  interferer  ASSR  =  30  dB 
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Figure  E-36.  Generic  signal  bounds  -  2  complex  Gaussian  signals.  100  looks.  SOI 
ASSR  =  50  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-37.  Const  ant -envelope  signal  bounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASSR  =  10  dB.  interferer  ASSR  =  10  dD. 
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Figure  F-.'IH.  Constant-envelope  signal  hounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASS'R  =  10  dB.  interferer  AS.XR  =  30  dB. 
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Figure  E-39.  Constant-envelope  signal  bounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASSR  =  10  dB.  interferer  ASXR  =  50  dB. 
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Figure  E-40.  Const  ant -envelope  signal  bounds  -  2  random  constant -envelope  signals. 
10  looks.  SOI  ASNR  =  30  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-41.  Constant-envelope  signal  bounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASSR  =  30  dB,  interferer  ASNR  =  30  dB. 
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Figure  E-42.  Constant -envelope  signal  hounds  -  2  random  constant-envelope  signals. 
10  looks.  .SOI  ASNR  ^  30  dB.  interferer  ASNR  =  50  dB. 


6f) 


116880-45  116880-44 


LOWER  BOUNDS  ON  DF 


SOURCE  SEPARATION  (BW) 


CE  SIGNAL  BOUND 
RANDOM  CE  SIGNALS 
UNIFORM  LINEAR  ARRAY 
10  ANTENNA  ELEMENTS 
2  SIGNALS 
10  LOOKS 
PI  =50dBASNR 
P2=  lOdB  ASNR 
NO  ARRAY  ERRORS 
-30  dB  ARRA  Y  ERRORS 
-25  dB  ARRA  Y  ERRORS 
-20  dB  ARRAY  ERRORS 
-15  dB  ARRA  Y  ERRORS 


Figure  E-43.  Constant -envelope  signal  bounds  -  2  random  const  ant -envelope  signals. 
10  looks.  SOI  ASNR  =  50  dB.  interferer  ASNR  —  10  dB. 
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Figure  E-44.  Constant -envelope  signal  bounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASNR  =  50  dB.  interferer  ASNR  =  30  dB. 
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Figure  E-45.  Constant-envelope  signal  bounds  -  2  random  constant-envelope  signals. 
10  looks.  SOI  ASyR  =  50  dB.  interferer  ASNR  ^  50  dB. 
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figure  E-46.  Constant-envelope  signal  bounds  -  2  random  constant -envelope  signals. 
100  looks.  SOI  ASNR  -  10  dB.  interferer  ASNR  =  10  dB. 
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Figure  E-47.  Const  an  t-erii'eJope  signal  bounds  -  2  random  constant-envelope  signals, 
lot)  looks.  SOI  .AS\R  =  10  dB.  interferer  ,4S.Vi?  =  30  dB. 
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Figure  E-48.  Const  ant-envelope  signal  hounds  -  2  random  constant -envelope  signals. 
100  looks.  SOI  ASXR  =  10  dB.  interferer  ASNR  =  50  dB. 
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Figure  E-49.  Constant -envelope  signal  bounds  ■  2  random  constant-envelope  signals. 
100  looks.  SOI  ASNR  =  30  dB,  interferer  ASNR  =  10  dB. 
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Figure  E-50.  Const  ant -envelope  signal  hounds  -  2  random  constant-envelope  signals, 
100  looks.  SOI  ASNR  =  30  dB.  interferer  ASNR  —  30  dB. 
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Figure  E-51.  Const  ant -envelope  signal  bounds  -  2  random  const  ant -envelope  signals. 
100  looks.  SOI  ASSR  ^  30  dB.  interferer  ASNR  50  dB. 
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Figure  F-52.  Const  ant -envelupr  signal  bounds  -  2  random  constant-envelope  signals. 
100  looks.  SOI  .ASNR  =  50  dB.  interferer  /  SNR  =  10  dB. 
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Figure  E-o3.  Constant-envelope  signal  hounds  -  2  random  constant-envelope  signals. 
100  looks.  SOI  ASSR  =  50  dB,  interferer  ASNR  =  30  dB. 
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P2  =  50  dB  ASNR 
NO  ARRAY  ERRORS 
-30 dB ARRAY  ERRORS 
^5  d^RRA  Y  ERRORS 
-20 dB ARRAY  ERRORS 
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Figure  E-54.  Constant -envelope  signal  bounds  -  2  random  constant-envelope  signals. 
100  hxtks.  SOI  ASNR  =  50  dB.  interferer  ASNR  =  50  dB. 
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